Critical behavior of loops and biconnected clusters on fractals of dimension d < 2. 
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, We solve the 0(n) model, defined in terms of self- and mutually avoiding loops coexisting with 

^ • voids, on a 3-simplex fractal lattice, using an exact real space renormalization group technique. As 

' the density of voids is decreased, the model shows a critical point, and for even lower densities of 

, voids, there is a dense phase showing power-law correlations, with critical exponents that depend on 

n, but are independent of density. At n = —2 on the dilute branch, a trivalent vertex defect acts as 
a marginal perturbation. We define a model of biconnected clusters which allows for a finite density 
of such vertices. As n is varied, we get a line of critical points of this generalized model, emanating 
(~| . from the point of marginality in the original loop model. We also study another perturbation of 

^ ' adding local bending rigidity to the loop model, and find that it does not affect the universality 

O ' class. 



Y . PACS numbers: 05.45.Df, 64.60.Ak, 05.50.+q, 64.60.Cn 

C/2 . I. INTRODUCTION 

The loop model is a very important model in statistical physics. It was defined originally in terms of the high 
temperature expansion of the n- vector model [H 0, The cases of n = 0, 1, 2 correspond to well-studied cases of 
■ self-avoiding polymers [3], the critical Ising and the XY models respectively. The model has been studied quite 
^ ' extensively in d = 2 dimensions, in several variants, including fully packed or dilute versions, and with loops of more 
O . than one type. One can determine the critical exponents of the model on the hexagonal and square lattices using 
the Bethe Ansatz technique, the Coulomb gas method, and numerical techniques involving exact diagonalization of 
transfer matrices of systems on finite width cylinders 0, [1, 0, 0, H, S • The critical behavior of the model is 
also related to the dimer model, edge coloring model, compact polymer model (for n < 2) [13, UH) > ^^'^ t;he hard 
K*" ' hexagon model (for n > 2) fisj. In the presence of a staggered field, the model gets related to the critical Potts model 

The loop model to be studied here is defined by the partition function 
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■^loop = ^n^W^, (1) 



where the summation is over configurations C of self- and mutually-avoiding loops on a lattice. In the above, n is the 
weight of a loop, lo is the weight of a vacancy (i.e., a lattice vertex not visited by any loop), C denotes the number of 
loops in a given configuration, and V the number of vacancies. 

In Eq. Il]), large values of uj correspond to a small average density of loops, and thus to the high-temperature phase 
of the n-vector model. As lu is decreased, the average density of loops and the mean loop size increases, and diverges 
' as w tends to an n-dependent critical value Lu^^(n). The critical behavior of the loop model at a;*+ gives the critical 
behavior of the n-vector model at its critical point in d > 2. For to < u;*+ and n < 2, the model shows a critical phase, 
which is called the dense phase. In d = 2, the critical exponents of the dense phase have been determined exactly 
[1, 0) B Hi • These vary with n, but are independent of the precise value of uj. 

The critical behavior of the loop model, and exponents of the dense phase have been less studied for d 2. One 
expects that the dense phase of the loop model for d > 2 would be related to the the low-temperature phase of the 
n-vector model. The latter shows power-law correlations in the entire low-temperature phase because of the existence 
of gapless Goldstone modes. In c? = 2, Jacobsen et. al. [l^ have shown that allowing loops to intersect or not leads 
to different critical behavior, and the dense phase of the loop model is different from the Goldstone phase. 

It therefore seems interesting to study the loop model in dimensions other than d = 2. In this paper, we study 
the loop model on fractal lattices with finite ramification index. We will take the 3-simplex lattice (see Fig. [T]) as the 
simplest example of this type. The treatment is easily extended to other fractals. The loop model shows a nontrivial 
critical behavior for n < 1 on this fractal, and we determine the critical behavior near the dilute critical point 1^*+. 
We also study the critical properties of the dense phase. Further we define a generalization of the loop model, to 
be called the biconnected clusters model — or "bicon clusters model" for short — in which we allow the summation in 
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FIG. 1: The construction of the 3-simplex lattice is shown from level r = to r = 1 and recursively for higher r. 
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FIG. 2: (a) A typical configuration of the loop model, (b) A typical graph showing several biconnected clusters that contribute 
to the partition function of the bicon clusters model. 

the partition function to include all biconnected clusters, as shown in Fig. [21 A cluster is called bi-connected, iff 
between any two points of the cluster, there are at least two disjoint self-avoiding paths. Clearly all simple loops are 
biconnected. It turns out that this changes the critical behavior of the model, and we determine the new critical 
exponents. 

A d — 2 variant of the bicon model is known as the "net model" (l6j . This model is relevant for the study of quantum 
models whose ground states are endowed with topological order. More generally, the constraint of /c-connectedness 
has also been imposed in classical models of clusters, such as percolation [Tt'I and the Potts model p^ . 

There have been several studies of critical behavior of statistical mechanical systems on finitely ramified fractals 
p^ . The Ising model was the earliest to be studied on fractals It was followed by a study of self-avoiding 

polymers on a 3-simplex lattice pill . Later self-interacting self-avoiding polymers [2^ , and other trails [23| , have been 
studied. The Lee Yang edge singularity for Ising model was considered in [23| . The collapse transition of branched 
polymers was studied in [2^. The distribution of sizes of erased loops for loop-erased random walks was studied in 
[2^ . In [27II , it was studied how the number of self avoiding rings going through a site varies with the position of the 
site. Polymers with bending energy were studied on 3-simplex and other fractal lattices in [28l |. For a recent review, 
see [29]. 

The detailed plan of the paper is as follows. In Section |TT] we find the critical branches of the loop model, and in 
section Hill the exponents are calculated. In section HVl we study the generalized "bicon model", as mentioned above. 
In section IVl we study the loop model with extra energies for local bending. In section IVTl we summarize our results. 



II. THE CRITICAL BRANCHES OF THE LOOP MODEL 



The recursive construction of the 3-simplex fractal lattice, through a series of levels r, is illustrated in Fig. [2 Let 
Zr be the partition function at level r, i.e. the sum of statistical weights of all configurations with loops contained 
within level r, and denoted by an empty triangle in Fig. |3l We may have loops which close at levels higher than 
r. Such a loop will come out of the corner vertices of an r-th level triangle. The restricted partition function for 
configurations with such an open chain passing out of the corner vertices of a r-th level triangle is denoted by Br and 
shown in Fig. [31 




Z B A C D 

FIG. 3: The schematic representations of various restricted partition functions Z, B, A, C and D for an r-th level triangle. 
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These functions Zj. and can be recursively related to their counterparts at the [r + l)-th level by the following 
equations: 

Zr+l = Z,^ + nBr^ (2) 

-Br+l = ZrBr^ + B^^ . (3) 

At the 0-th level the lattice has just a single site. In that case Z^ — lo is the weight of it being empty, and i?o = 1 is 
the weight that a loop passes through it. If we define ujr =^ Zj./ Br, then the recursion relation for the latter is 

Lol + n 

^r+l = . (4) 

ijJr + 1 

At the 0-th level, u)q = Z^j B^ = u). The total free energy on an infinite lattice is 

F{n,uj) = lim ^ In Zr+i. (5) 

r — >oc ^5 

Eq. ^ above may be rewritten, using Zr+i = Z'^{Zr+i/Zi^), and after a little algebra we obtain 

F{n,u;)^-ln{l + nco-^) + -Fin,—^j, (6) 

which is the usual form [30| of the recursion relation for the free energy under the real-space renormalization group 
(RG). Apart from = oo, the finite- valued fixed points of the RG fiow ^ are the roots of the following equation: 

Lul-ujl-Lj,+n = 0. (7) 

The special case n = corresponding to self-avoiding polymers was studied earlier [2l| . 

Eq. ([7]) has three real solutions for —5/27 < n < 1, and only one real solution for the regions n > 1 and n < —5/27 
(see Fig. U]). The three fixed points for the region —5/27 < n < 1 are given by the formulae 

n = JL(ll-16cos(0)), (9) 

where 6 e [0, Bvr]. The three solutions correspond to the following subdivision of the parameter range of 9: 

• For 9 e [0,7r], n runs from —5/27 to 1, and uj from 5/3 to 1. We shall refer to this line of fixed points as the 
dilute branch and denote it by Lu^-i- (n). It is physical (n > 0) for 9 > cos-i(ll/16) ~ 0.8128. 

• For 9 G [tt, 27r], n runs from 1 to —5/27, and lo from 1 to —1/3. We shall refer to this as the dense branch and 
denote it by It is physical (n > and w > 0) for < 27r - cos"^(ll/16) ~ 5.4704. 

• For 9 e [2tt, 37r], n runs from —5/27 to 1, and oj from —1/3 to —1. 

This unphysical critical point at negative uj may be called the Yang Lee edge singularity in our problem. We 
shall denote it by ujyl- However, we note that the critical behavior at this "Yang Lee edge singularity" is not 
independent of n. 

The three branches a;*_ and luyl are shown in figure [4] Note that unlike the two-dimensional loop model, 
where the critical region between the dilute and dense branches extends up to n = 2, the domain of criticality of the 
present model is smaller. But here too, the upper branch is a repulsive fixed line. Starting from any lu > a;*+ the RG 
flows take one to lu = oo while starting from any uj < the lower branch a'*_(n) serves as a line of attraction of 
the RG fiow. We will refer to the region with w < ti;,+ as the dense phase. 

For n < —5/27, the dilute branch can be extended backwards by replacing 9 = i9 (with 9 E [0,oo)), i.e replacing 
"cos" by "cosh" and 9 hy 9 in Eqs. ([5]) and Q. Although this regime is unphysical, we will show in section UlI DI that 
the point {n = —2, w = 2) on it is of some signficance. 

For n > 1, the nontrivial fixed point other than tj, = oo is negative, and is obtained by replacing 9 — + i9 (with 
6 E [0,oo)), i.e. replacing "cos" by "— cosh" and by 6* in Eqs. ^ and There is no nontrivial critical point for 
positive u). It can be shown [31^ that the loop model for n — 2 can be exactly mapped to a weighted 6- vertex model, 
which also can have no nontrivial critical point for positive weights in this regime. 
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FIG. 4: A schematic picture of the three hnes of fixed point 1.0^+, uj,- and ujyl plotted against n. 



III. THE EXPONENTS OF THE LOOP MODEL 



In two dimensions, for a generic point lying on the dilute branch or in the dense phase, the loop model is known 
to have an infinite set of exponents characterizing its critical behaviour [l^ . Is^ . The full set of principal exponents 
consists of an infinite number of correlation exponents associated with defects called the fc-string defects, plus one 
thermal exponent. 

A /c-string defect has an integer k(> 1) number of "open chains" originating from a local neighborhood and 
terminating in an anti-defect very far away. Since chains are disallowed in the original loop model, such structures 
may be imagined to be created in isolation, as defects, upon application of suitable small external fields. A closed 
loop with two points marked on it is equivalent to a 2-string defect, and this observation can be used p^. [32^ to 
find the fractal dimension of the closed loops in the loop model. Although there is no upper bound on k for ordinary 
Euclidean lattices, on fractals due to finite ramification, the exponent spectrum becomes finite. For the present model 
on the 3-simplex fractal lattice we can have such defects only with k < 3. [On the Sierpinski gasket, one could have 



Apart from the fc-string defects, the loop model also has the thermal defect and different exponents associated with 
it, namely the specific heat exponent, the correlation length exponent, etc. Among the latter only one is independent 
and others can be related to it by identities. Thus there are four independent exponents for the loop model on 
the 3-simplex lattice, three /c-string exponents plus one thermal exponent, and we calculate them in the following 
subsections ImHImP] 



Let us define (uj — lu^j^) as the "thermal distance" from the critical point lj^,^. Starting from luq = tu > co^^, under 
succesive iteration of Eq. the separation (5^ = {LUr—'-^*+) increases with increasing r. We note that (5o — {w — uj^^). 
Using Eqs. ^ and ([7]), we easily find that Sr = X^^o with 



The thermal exponent yx is defined as usual by the fact that at the r-th step of RG transformation, the thermal 
distance Sr ~ (2^'^)''(5o. Thus comparing the two expressions for Sr we obtain 



k < 4.] 



A. The exponents related to thermal excitation 




(10) 
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The specific heat exponent a may be defined via the scaling behavior of the non-analytic part of the free energy, 
namely F {u — uj^,-^)'^~°' . Substituting the latter scaling form into Eq. ^ and comparing the singular terms, a 
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straightforward calculation yields 

In 3 



In At ' 



(12) 



where At is given by Eq. ([TC 

Thermal fluctuations create thermal defects and the correlation function C(R) for two such defects separated by 
distance |R| is defined for u > as 

The correlation length ^ is finite for uj > ijJ*+ and the scaling relation ^ ~ (a; — uj^^)~'^ defines the correlation length 
exponent v. Under RG, if ^ ~ 6, since (uj — a;*_|_) ^ b^^^ it follows immediately that 

^ = i/yT-T^. (14) 

In At 

For to — in Eq. (|13p . ^ — > oo and c = 2xt, where xt is called the thermal correlation exponent at criticality. 
From standard RG arguments [s^l it can be shown that xt — d — yr, where 

In 3 , , 

is the box dimension of the 3-simplex fractal lattice. From Eqs. (fT^ and we see that the hyperscaling relation 
dv = 2 — a holds. 

Apart from the exponents of thermal origin discussed above, there are several other exponents caracterizing the 
behavior of other physically interesting observables, which can be expressed in terms of the exponents described above: 

• The fractal dimension of the loops given hy df relates the length s and radius i? of a loop as s ^ W^f . The 
exponent d/ is distinct for the two critical regions uj = a;*+ (dilute branch) and uj < (dense phase). We 
show in section IIII CI that they can be easily related to the 2-string correlation exponent X2 ■ 

• For Lu > the distribution of loop sizes s has an exponential cutoff ~ exp(— as). The cutoff length a^^ ^ (^"^f , 
where f is the finite correlation length in the disordered phase. One can find a^^ as in the derivation of v above. 

• For the critical regimes lo — lu^^ and lo < lo*+, the probability distribution of loop size s is an unbounded power 
law ~ s~^'^~^\ The exponent r is distinct for uj — uj^,+ and uj < uj^^, and it is easy to express the exponent r 
in terms of the exponents X2 and the fractal dimension d/ of the loops (see section IIlIC|) . 

The exponents like df and r can all be expressed in terms of the "fc-string defect" exponents which we derive in the 
following three subsections. 

B. A 1-string defect 

Let the application of a small external magnetic field hi create an open chain, or a 1-string defect, with the 
magnetic scaling exponent yi associated with it. Assuming for hi — the system is critical, a finite hi introduces a 
finite correlation length ^ given by 

ahi)^h;'/y\ (16) 

Eq. ([TB| defines the scaling exponent yi. In real space RG, upon coarse-graining from level r to level (r -I- 1), ^ 
increases by a factor of 2, while hi increases by a factor of Ai, i.e. ^{hi) — 2^(Ai/ii). Below, we will calculate Ai, and 
using this with Eq. (|16p. one gets 

The probability G(X — Y) [s^l of k strings originating at X and ending at Y in the critical phase is given by 

G'(X- Y) - 1/|X- Yp^". (18) 
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^From standard RG analysis it is known [3u| that the correlation exponent Xk is related to yk via the relation: 

Xk^d-yk- (19) 
Thus for our 1-string defect, by using Eq. (fT7|) and box dimension (|15p. we get 

ln(3/Ai) 



Xi 



yi 



In 2 



(20) 



Now we proceed to obtain the scale factor Ai. First we note that an 1-string defect may have its one endpoint 
inside an r-th level triangle. The corresponding restricted partition functions Ar and Cr are represented by diagrams 
A and C in Fig.[3l The recursion relations for Ar — Ar/Zr, Cr = Cr/Zr, and Dr = Dr/Zr are as follows: 



Ir+l 



Aril + 2lu-^ + 2cj-2) + CrLO-^in + 2) 



1 + nUJr 



Cr+l — 



{Ar + 'iCr)0J'r 



1 + mO' 



-3 



D 



r+l 



Al + Cluj-^{n + 6) + AArUj-^Cr + 2Aluj-^ + br{2uj-^ + Zoj-^) 

1 + ntOr'^ 



(21) 
(22) 
(23) 



We can linearize Eqs. (PT|) - ((^ around the fixed point {A* ,C* ,D*) = (0, 0, 0), and we are left with the matrix 

a + 2 a(n + 2)6^^2 

-3 



(24) 



with a — 1/(1 + nuj^ ). The largest eigenvalue of the matrix is Ai and is given by 

Al = ^ (Sa + Lol{a + 2) + Va2(17 + An) - Qauul{a + 2) + ijjt{a + 2)^) . 



(25) 



Note that in Eq. (j25p . must be replaced by uj^j^ or in order to describe the dilute or dense branch, respectively. 
Using Eq. ((25l) . we may thus write Ar ~ A^Aq. It is clear that is proportional to the small external field hi, and 
thus using Eq. 125]) . we may read off the values of yi and xi in Eqs. pT]) and ([20|l . 

Another exponent of interest is 71 associated with the approach from above of the dilute critical branch. It is 
defined via the scaling behavior of the average open chain length (Zi) ^ [uj — In general. 



(26) 



where c„ is the number of distinct configurations each with an open chain length n in a lattice of size iV, and Pn is 
the relative weight factor for such a configuration normalized by the partition function. On our fractal lattice, this 
becomes 



I {?>Al_i + ■iuo~\Al_i + ■iuj:;l^br-i 



E 



^r 



1 + nuj'r 



(27) 



where the statistical weight of fully containing an open chain at the (r + l)th level is fr+i — 'iA^Zr + 'iBrA^. + 3BrDr- 
For a fixed and small distance from the critical point So — (t-^ ^ ^*+)) since (h) is finite, the sum in Eq. (|27p is sharply 
cut off at some finite level r = rg. We note that for r > tq, Ai ^ ^ , and Dr ^ X^^° , and ^0 ^ (constant)/A^'' (see 
Eq. (fTO|)'). It immediately follows from the scaling relation (li) ~ ''^ that 



71 



ln(AV3) 
In At 



As a check of consistency, we may verify the following exponent equality |32l |. using Eqs. 

ji — [d — 2xi)v. 



, (EOD and H 



(28) 



(29) 
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C. A 2-string defect 



Let a triangle at the r-th level which has endpoints of 2-strings on two neighboring sites inside it, be defined to 
have a statistical weight D^. Note that this weight is different fi'om Dr (see section UlI CI above), since the latter puts 
no restriction on the location of the endpoints. Further defining D = D / Z , we easily see that its recursion equation 
is 

1 + nUJr 

Replacing Ur with the fixed point value in the above equation, we get D^. ^ ^2^0-1 with 

2w,-i + 3w,"2 + 3 

- ^3 - --Ti l-^lj 

1 + ncj* cj, + 1 

where we have used Eq. ([7]). Note again that in Eq. (j3ip for A2 we have to use w» ~ uj^^ for the upper branch, and 
Lu^ — for the lower branch. If a small field creates a 2-string defect and the finite correlation length arising 
due to it is ^(/i2) ~ /i^^^^^, we find that the scaling exponent 1/2 for a 2-string defect is 

2/2 = lnA2/ln2. (32) 

^From Eq. (|19p we conclude that the corresponding correlation exponent X2 is given by 

X2^d-y2. (33) 

It is important to note that 1/2 is also the fractal dimension df of the loops. This is explicitly seen as follows. Let Ir 
be the typical length of a segment of a loop that goes from one corner vertex to another in one i3-type triangle (see 
Fig. [3]) of order r. Then the typical length of a loop which closes at the (r + l)-th level is Sr — Sir- The recursion for 
Ir is: 

The above Eq. follows from the fact that at level (r+l), a typical loop length Ir+i can be made of two Ir segments 
with statistical weight ZrB^/Br+i, or three Ir segments with statistical weight B^/Br+i 33]. On the other hand, the 
typical diameter Rr of such a segment of length Ir is given by i?^ ~ 2*^ . Combining these two results and using the 
definition of d/ given by the scaling behavior ~ ~ Rr , we conclude that 

We note that df, through A2, is distinct for — (jJ*+ (dilute loops) and — cj*- (dense loops). For n — and 
in the dense loop phase, A2 — 3, and so df coincides with the box dimension d of Eq. (llSp . This was to be expected, 
since the limit of dense loops (w^,,- 0) in fact means that there is a single loop covering the entire lattice, i.e., the 
loop is Hamiltonian. For dilute loops, we have ijJ*+ ~ ^(1 + a/S), whence A2 = ^(7 — and df ~ 1.2522. 

The probability distribution P{s) of loop size s is an unbounded power law ~ s~('^~^) for both the dilute branch 



and the dense phase. The exponent r can be related to X2 and df following a general derivation as in [3j|. Let G's(R) 



be the probability that two points separated by R are on a loop of size s. The expected scaling form of G's(R) is: 

G.(R) ^ .-|Rr-/r . (36) 

Firstly, the sum of Gs(R) over all the space, i.e. / d'^RG's(R), is nothing but number of points on the loop = s. This 
gives a relation d — ci — df{l — m). Secondly, by definition (see Eq. (HH])) we have G(R) = / dsP{s)Gs{'R.) ^ jRj^^^^. 
The latter relation combined with the former gives 

df{3-T)=d-2x2, (37) 



which is the desired result. 
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D. A 3-string defect: 



Let a triangle at the r-th level which has endpoints of 3-strings on three neighbouring sites inside it, be defined to 
have a statistical weight E^.. The recursion relation of = E /Z is 

3^-2 

£^,.+1 = , _3 E^. (38) 

1 + nUJr 

Assuming that a small field creates a 3-string defect and the corresponding finite correlation length ^(/is) ^ '^Iv^ ^ 
we find that 

2/3 = lnA3/ln2, (39) 

where 

A3 = (40) 

The correlation exponent 2:3 corresponding to 2/3 is given by Eq. (|19p with fc = 3. 

We now note something very interesting. If we extend the upper dilute critical branch to negative values of 
n < —b/27, we find that at the special point (w, = 2, n = —2), we have A3 = 1 and 2/3 = (see eqs. and PO)) ). 
Thus although the 1-, 2- and 3-string defects are relevant for general n on both the critical branches, at the point 
(w, = 2,n = —2) the 3-string defect becomes marginal. In the next section HVl we define a more general model called 
the "biconnected cluster model" which allows for vertices of degree 3, each occuring with a finite weight ^3. We find 
that the latter model has a new critical line in its larger parameter space, which precisely meets with the line of 
critical points corresponding to = at the point (w* = 2, n = —2) on the extended dilute branch. 



IV. THE BICONNECTED CLUSTER MODEL 



The finding in section [III Dl that 2/3 = at a point in the (w,?t.) parameter space, motivates our study of a model 
called the biconnected cluster model in which the allowed vertex degrees are 0, 2 and 3. In other words, 3-strings are 
allowed to emerge from any vertex (see Fig.lSfa)). The connected components are further required to be 2-connected, 
i.e., they cannot be disconnected upon cutting a single link (see Fig. [Sj^b)) — note that this is a stronger requirement 
than simply disallowing vertices of degree 1. Henceforth we will refer to this model as the "bicon clusters model", 
to distinguish it from the "loop model" studied this far in the paper. A similar model — with no requirement of 
biconnectedness — was studied in two dimensions in ^16] under the name of "net model" . The partition sum for the 
bicon model is 

C 

where the summation is over configurations C of any number of self-, and mutually-avoiding biconnected clusters. 
Here n is the weight of a cluster of any size, u> is the weight of an empty vertex, and uj^ is the weight of a vertex 
of degree 3. Further, C denotes the number of clusters, V the number of vacancies, and U the number of vertices of 
degree 3, in a given configuration. Note that by setting u)^ = 0, the bicon clustes model reduces to the usual loop 
model. Some of the configurations that we exclude from the model (see Fig. [5l^b) for an example) will be treated as 
defect configurations called "fc-defects" in section HV Bl below. 

The three possible vertex configurations with weights w, 1 and ujj, are shown in Fig. [51 Note that if Mk is the number 
of vertices of degree k, we have the simple topological identity 2A2 + iM^ — 2L, where L is the total number of links 
in the configuration. Therefore As ^ U is necessarily even. Accordingly, the bicon clusters model partition function 
is a function of (ws)^. At the r-th level, the real space RG closes for three restricted partition functions schematically 
shown as Z, F and G in Fig. [S] While Z^. is the partition function summing over configurations with no strings 
coming out of the corner vertices, Fr and Gr have two and three strings coming out of the corners, respectively. Note 
that the constraint of biconnectedness implies that it is not possible to have a configuration with one string coming 
out of the corner vertices. At the level r = 0, 



Zq — uj; Eq — 1; Go — W3. 



(42) 
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FIG. 6: The three vertex weights for the bicon clusters model are shown in the top row. In the bottom row, the three schematic 
representations of the restricted partition functions are shown. 



At the r-th level, the recursion relations are: 

Zr+i = Zl + nF^ (43) 

Fr+l = +F^Zr + FrGl (44) 

Gr+i = Gl + 3GrF^ (45) 
Further defining Fr = Fr/Zr, and G,. = Gr/Zr, we get 

^ 1 + nF? ^''^ 
~ G\ + 3GrF^ 

Note that Fq — uj~^ and Gq — lu^u;^^. There are several fixed points (i^*,^?*) of Eqs. (|lS)) - ([Tf|) . Apart from the 
trivial weak coupling fixed point (0, 0) and the fixed point (cxd, oo) corresponding to all bonds being fully covered, 
there are the dilute and dense fixed points of the loop model which we denote by (F+ , 0) and (F_ , 0) , respectively. 
But most importantly there is a nontrivial fixed point (F^ , G, ) given by, 



F» = 1/2 and G, = ^ ^ + ^. (48) 

The latter defines a new critical line in the (w, 71,^3) space, which terminates on one end at the point (lu — 2,n = 
—2,uj3 = 0); note that this is the point where we found ys to be marginal for the loop model in section [III D I The 
dense-phase fixed point (i^_, 0) is unstable to introduction of trivalent sites, and the limiting behaviour of the critical 
net model is governed by the fixed point (F^, G*). 

If the starting value of Fr and Gr is near F^ and G», but a bit larger, it is easy to check that Fr and Gr diverge 
to infinity as r and r^/^ respectively. This implies that in the dense phase of the bicon clusters models, corner sites 
of triangles of high order belong to the infinite cluster with a large probability, and this probability tends to 1 as r 
tends to infinity. 



1- defect 




2- defect 
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FIG. 7: The 1-, 2- and 3-defects are shown along with their shorthand representations (filled black blobs). The restricted 
partition functions contributing to the three defects are shown schematically on the right. 



A. The fractal dimension of the biconnnected clusters 



Just like in the loop model, we would like to find the fractal dimension of the clusters in the bicon clusters model 
at its critical point. Consider two marked points A and _B on a cluster, free of tadpole overhangs. Imagine that these 
marks are "defects" created by some external field h. In addition to F and G, we now define new functions F"^ and 
G™ which are analogous to F and G, except that they represent configurations in which there is exactly one marked 
point. At level r = we have: 

F^ = h 

= huj3. (49) 

The recursion relations are as follows: 

= 2FrF;''Zr + 3F^F^ + 2FrGrG'^ + GlF^ (50) 
G*"!^]^ — QF^F^'^Gr 4" 3F^G^ + 3G^G^ . (^1) 

The largest eigenvalue of the 2x2 matrix obtained by linearizing the above equations around the fixed point of 
Eq. gll) is 

^bicon ^ 2(7 + 7i) + Vn2 + 20n + 52 ^ 

and the fractal dimension of clusters on the new critical line is 

d^'™" = ln(A^'™")/ln2. (53) 

The scale factor Aa"^"" A2 = 7/3 and thus d^'™" ^ df, the fractal dimension of the loops (see Eq. (|35p ). at the 
point (a; = 2, n = — 2, 073 = 0), as expected. 



B. The fc-defects 



What kind of defects are natural extensions of fc-string defects, appropriate to the net model? A possibility is what 
we call the "fc-defects" shown in Fig. [T] A 1-defect has multiple clusters connected in a series by strings, and the two 
dangling ends are marked; they are depicted in short-hand as filled black blobs (see Fig. [7]) . The obvious motivation 
for defining such defects is that as W3 ^ 0, they become our usual "fc-string" defects in the loop model. The restricted 
partition functions contributing to each fc-defect is shown alongside the defects in Fig. [Tj 
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n 


-2 


-1 





1 


2 


oo 


^bicon 


10/3 


3.23 


3.21 


3.16 


3.13 


3 



TABLE I: Some values of Ai"^°" versus n associated with the 1-defect. 



1- defect dimension: In Fig. [71 A ,C ^ Xi and X2 represent the relevant restricted partition functions for the real 
space RG of a system containing a 1-defect. At the 0-th level, A^^ = h (where h is an external field which gives rise 
to this defect) and Cq = ATi o = X2,a = 0. The recursion equations for A , C , A"i and X2 (scaled by Z as usual) are 
as follows: 

r' _ A',{l + 2Fr + 2F^) + {2 + n)F^Cl + 2GrFrXir+F^X2r 

'■+1 ^ T"^ £^3 ^ ^ 

1 + nF^ 

= ^^^^ 

~ 2GrFrX^ + AGrFrCl + {2Fr + SF^ + (^Z)^^ ^ _^ 2GrFrX2,r 

= (56) 

- 6GjCl+6GrFrXi,r + {SF^ +3Gl)X2,r 

^^-+1 y:^. 

Solving the 4x4 matrix obtained by linearizing the above equations about the fixed points in Eq. (|48p . the largest 
eigenvalue A^"^"", gives the scaling dimension yY'^°'^ associated with a 1-defect: 

^bicon ^ A'l-'^'Vln 2. (58) 

Some values of A^"^°" versus n are given in Table 1. Note that upon approaching the point [n — — 2,cj — 2,ujj, — 0), 
we have ytucon _^ j recover the scaling exponent of the 1-string defect in the loop model, as expected. 

2- defect dimension: Again the relevant restricted partition functions are shown in Fig.[7]by representative symbols 
X3 and X4. The recursions for the latter scaled by Z arc as follows: 

~ X,A2Fr + 3F^) 

iJr ^X^^rFj-Gr + (3F^ + Gl)Xir 

Xi^r+1 = ■ ; (60) 

1 + nF^ 

Assuming Xg^^ - {\f''°'^Y the scale factor A^'™'^ = {2uj-^ + iuj-'^)/{l + 7iuj-^) and independent of ^3- Putting the 
critical point value uo^ = 2 (Eq. ((48|) ). we get: 

A^'™" = (61) 
8 -I- n 

3- defect dimension: For this defect the relevant restricted partition function is shown in Fig. [7] and is represented 
by the symbol X^. The recursion relation for X5 = X^/Z is: 

X5,r+1 = 1^ (62) 

1 -I- nl'^^ 

and the corresponding scale factor on the critical line (Eq. 

;^bicon ^ _J_ /ggN 

3 n+8 ^ ' 



Thus for n > —2, the scale factor A3"^°" < 1 and hence the 3-defect is irrelevant on the entire new critical line 
(Eq. (O). This was certainly to be expected, since it is exactly this 3-defect that induces the flow from the 0^3 = 
line to the new critical line. At the point (71 = —2,uj = 2,uj^ = 0), we have Aj^"^"'^ = 1 as expected, since A3 = 1 for 
the loop model at that point. 
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Z„ = CO Bo= Xl Co = 1 Do= 
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FIG. 8: The different local configurations for a vertex is shown. For a covered vertex the loop can either go straight or bend 
in two ways. 




FIG. 9: Possible restricted partition functions are represented schematically above. Note that the external legs covered by 
loops are distinct in all of them. 

V. THE LOOP MODEL WITH LOCAL BENDING ENERGY 

The problem of self-avoiding polymers with bending energy has been of long standing interest, and was first 
introduced by Flory [35| . With high energy cost for bending, the polymer is in an ordered state (with minimal 
bending), while reducing the energy cost leads to a disordered (but critical) state. The nature of the phase transition 
separating the two phases was unclear for a long time, and finally it was shown recentl y th at for a compact (i.e., space 
filling) polymer on a two-dimensional square lattice the transition is continuous [s^. l37j|. The latter works actually 
dealt with the full loop model and obtained the relevant results for the polymer by taking the n — > limit. In a 
similar spirit we would look at the loop model for general n with local bending energy on 3-simplex fractal lattice; 
unlike (Sa. [37| our loops are not compact. We note that the n ^ limit has already been studied earlier in [11] where 
it was found that the bending energy is irrelevant and no new fixed points appear in the extended phase space. We 
show below that the same is true for general n. 

The model to be studied here is defined as 

^loop = ^n''a;''A]'^A^^ (64) 
c 

On a 3-simplex fractal lattice, at any vertex, a loop can go straight, or bend by or ^ (see Fig. [5|). Accordingly 
we define local vertex weights 1, Ai, A2 respectively, for the three cases. Correspondingly, the number of vertices 
covered with loops having ^ and ^ bends are respectively Vi and V2- The remainder of the symbols in Eq. ([M]) 
are as in Eq. ^ for the partition function of the loop model. For doing the real space RG in this case, we require 
four restricted partition functions with distinct external legs covered by loops as shown in Fig. [51 These are further 
represented by symbols Z,B,C and D, and at level 0, the values of these partition functions are shown in Fig. [S] 
Note that for Ai = A2 = 1, the diagrams B, C and D are equal. 
The recursion relations for the restricted partition functions are 



?r+l 

l2 



Cr+1 — BpDj, + ZrCj. 
Dr+l = BrCrDr + ZrCrDr, (65) 

with Zq = uj, Bq = \i, Co = I and Dq = X2 (see Fig. If we scale B, C and D by Z, we get 

~ Br+l _ BrC^ + D'^ 

Hr+l - ^ - -3 (66) 

1 + nB^ 

~ Cr+1 BrDl + Cr , . 

<-r+l - ^ = — — (67) 

Zr+1 1 -I- nB;^ 

~ Dr+l BrCrDr + DrCr , . 

Dr+l = ^ ~ (68) 

^r+1 1 -|- nBi: 



13 




FIG. 10: The RG flow shown on a plane containing the (1, 1, 1) direction in the B, C, D space for a fixed n. 

From the above Eqs. (pS)) -l|68 p we get the following fixed points (assuming Ai and A2 7^ 0): 

B^^C,^ and 1 + nC^ = + . (69) 

Comparing this with Eq. ([7|) we see that non-zero Ai and A2 have no new effect on the loop model and no new critical 
points come into being. The irrelevance of bending energy is shown in gray in Fig. 1101 The figure shows two fixed 
points in the (1, 1, 1) direction in the B,C,D space for fixed n. The attractive fixed point and the repulsive fixed 
point in the (1, 1, 1) direction are exactly the same as two points for a fixed n (belonging to the two branches) in 
Fig. m The RG flow diagram on a plane containing the (1,1,1) line is shown in Fig. [TOl 

The irrelevance of bending rigidity is not found for all fractal lattices. In [28j, as well as in our calculation above, 
exponents do not get affected on the 3-simplex lattice. Similarly, in (2^ it was shown that there is no effect of 
the strength of self-interaction on the swelling exponent of a self-avoiding walk on the Sierpinski gasket. However 
on another fractal, namely the branching Koch curve (BKC) the exponents do change. This was shown for the 
self-avoiding polymers [1^ and recently the full loop model for general n [3l| on BKC. 

VI. CONCLUSION 

In summary, we have studied the n- vector model on a fractal lattice with dimension d < 2. We have shown 
that — ^just like its counterpart in d = 2 — it has two (physical) critical branches, referred to as 'dilute' and 'dense' 
loops. However, while these branches exist for 0<n<2in(i = 2, on the 3-simplex lattice they are constrained to 
< rt < 1. We have explicitly characterized the critical behavior of this model, in terms of three fc-string and one 
thermal exponent. 

We note that the dense phase exists on this fractal for n < 1, and the dense phase exponents have a nontrivial 
dependence on n. This is in contrast to the expected behavior of the n-vector model for d > 2, n > 1. For n > 1 and 
d > 2, the low temperature phase of the n-vector model in zero field is characterized by a spontaneous magnetization. 
The entire low-temperaure phase is critical, with infinite correlation length, because of the presence of Goldstone modes 
in the system. There are n—1 Goldstone modes, as there are (n — 1) directions orthogonal to the magnetization. For 
each of these modes, the excitation energy for an excitation of wave number k varies as fc^. Using this, and the fact 
that the spin wave approximation becomes asymptotically exact at low temperatures, one can easily deduce results 
like the spontatenous magnetization at low temperature T varies as 1 — A{n — 1)T'^/^ for general n. Also, the mean 
energy at temperature T varies as E{T) = E{T = 0) -I- B{n — l)T''/^+^. In all these cases, only the amplitudes are 
proportional to (n — 1), and the critical exponents do not depend on n. However, for the fractal we studied, the 
dense phase exists only for n < 1. In this case, if we do a naive analytic continuation in n, the number of Goldstone 
modes becomes negative. Then, for n < 1, magnetization density becomes > 1, and E{T) decreases with increasing 
T. Clearly the argument that the exponents of the phase are determined by Goldstone modes no longer applies. 
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We have further generahzed the loop model into the 'biconnected clusters model' which allows for a finite density 
of 3-valent vertices. This model was found to support a new line of n-dependent critical points, emanating from the 
n = —2 point on the dilute branch in the original model. At this point, the 3-string defect is a marginal perturbation, 
which is clearly a necessary requirement for the emergence of a new critical line in the generalized model. 

An interesting line of future research would be to study models allowing for 3-valent vertices in two dimensions. 
Also in d = 2, the 3-string defect is marginal at the n = —2 point on the dilute branch. Within conformal field theory, 
the latter point is a theory of symplectic fermions with central charge c = — 2. One might speculate that in this case 
the perturbation could be exactly marginal, and hence generate a line of c = — 2 theories with continuously varying 
critical exponents. Such theories were recently shown to exist "38^, and can be produced [38| by perturbing the c — —2 
point on the dense branch (i.e., n = 0) by a finite density of sia;-valent vertices. 
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